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Abstract: N=4 Super-symmetric Yang-Mills theory (N=4 SYM) in the strong coupling regime has been 
^ Successfully applied (through the AdS/CFT correspondence) to the description of strongly coupled plasma 
which is a multiparticle state. Yet, the high-energy scattering in the strong coupling limit of N=4 SYM 
^ is purely elastic, so this multiparticle final state can never be produced: this is because in this limit the 
m ^ theory is dual to weak supergravity, and the dominant interaction is the elastic graviton exchange. Here we 
^ propose a resolution of this dilemma by considering the contribution of D-instantons in AdS§ bulk space 
to the scattering amplitude. We argue that D-instantons coupled to dilatons and axions are responsible 
for multiparticle production in strongly coupled N=4 SYM, and the corresponding cross section increases 
with energy. We evaluate the intercept and the slope of the corresponding Pomeron trajectory in terms 
of the typical size of the D-instanton, and argue that the resulting physical picture may resemble the real 
world. 
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1. Introduction 

N=4 Super-symmetric Yang-Mills (SYM) theory does not possess the properties of confinement and asymp- 
totic freedom. Therefore, it cannot be considered a theory adequately describing the strong interactions 
in the real world. However, this theory has two unique attractive features. First, N=4 SYM at weak 
coupling has some properties similar to perturbative QCD (see Refs. Q §]): the Operator Product Expan- 
sion (OPE) and the linear evolution equations for Deep-Inelastic Scattering (DIS) are valid, as well as the 
BFKL equation for the high energy scattering amplitude. The high energy amplitude reaches the unitarity 
limit, or the "black disc" regime in which half of the cross section stems from the elastic scattering and half 
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from the inelastic processes of multiparticle production. Second, in the strong coupling regime N=4 SYM 
can be reduced to weak classical gravity in the AdS$ bulk space through the AdS/CFT correspondence 
P], and can be solved analytically. Therefore, for physical systems for which neither confinement nor run- 
ning coupling are essential, the N=4 SYM theory can provide guidance in the strong coupling regime. In 
particular, it may lead to a better understanding of phenomena in strongly coupled QCD matter. Indeed, 
some of the main successes of N=4 SYM to date have been achieved in the description of the quark-qluon 
plasma at strong coupling || ||, 0] . 

dilate + axions 2 However, the physical picture of the high-energy 

collision in the strong coupling regime of N=4 SYM 
turns out to be completely different (8|, |9|, [IT], [T^, 
14] from the one at weak coupling: the multi- 



unpartides 



dilatons + axions 




particle production in this case is apparently unim- 
portant and the main contribution stems from elas- 
tic and quasi-elastic (diffractive) processes when the 
target (proton or a nucleus) either remains intact, or 
is slightly excited. This picture contradicts both the 
QCD expectations |l|, [E§, [L7|, @ |l|, 
able experimental data. 



and the avail- 



Figure 1: The process of multiparticle production in 
N=4 SYM 

Since the creation of quark-gluon plasma neces- 
sarily proceeds through the multiparticle production, we face a rather controversial situation. Indeed, the 
multiparticle systems which we hope are adequately described by the strongly coupled N=4 SYM are not 
produced in this theory. The reason for this is simple: the main contribution to the scattering at high 
energies stems from the exchange of the graviton which has the spin equal to 2 and generates large but 
real scattering amplitude. The unitarity constraint for the scattering amplitude has the form 

2 



2lmA(s,b) 



\A(s,b) \ 2 + O 



1.1 



where A = gy M N c = 4irN c g s . We consider the regime where the number of colors is large N c S> 1, 
the SYM coupling is strong A > 1, but the string coupling constant is small j s < 1. This allows for a 
perturbative dual description in terms of the graviton interactions. Geometrically, the radius of the AdS§ 
space R is related to the coupling constant by R 2 = vA a' (2ira' is the inverse of the string tension) so 
that large values of A 3> 1 correspond to the AdS$ space of large radius and small curvature for which the 
perturbative gravitational description is valid. 

The unitarity constraint of Eq. ( |1.1| ) follows from the eikonal formulae (see Refs.|| |l^, 12, [14|]) that 
sum the main contributions of the order of (g s a' s) n and describe the scattering amplitude in the kinematic 
region where g s a' s ~ 1. The graviton production in the central rapidity region leads to the cross section 
which is proportional to g 6 a' 2 s 2 (see Fig. §-a ) and is much smaller than the eikonal diagram of Fig. §-b 
describing elastic scattering. The estimates for the smallness of the multiparticle production cross section in 
Eq. fll.l| ) follow from the fact that the reggeized graviton has a trajectory with the intercept jo = 2 — 2/y/X 
(seeRef. §). 
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Figure 2: The graviton production in N=4 SYM. Gravitons are denoted by wavy lines while dilatons by dashed 
ones. 

Therefore at g s -C 1 the emission of gravitons is suppressed and, at first sight, we are doomed to a small 



multiparticle production cross section that follows from Eq. (1.1). However, this conclusion is premature 
since the production of a large number of particles (n 3> 1) may not be adequately described by the lowest 
order perturbation theory. Instead in this case we can rely on the classical (semi-classical) approach. 
Indeed, the phase-number uncertainty principle An Acfi > 1 implies that the variance of the phase <p of the 
field describing the large average number n of radiated quanta should be small: Acft 2 oc 1/n <C 1, so the 
field is semi-classical. 

The deep connection between the classical solutions and the multi-particle production has been realized 



long time ago pl| , 22, [24[. Let us briefly remind the main ideas underlying this connection. Consider 
first a perturbative expansion of some amplitude F in powers of the coupling constant e: F(e 2 ) = ao + 
a2e 2 + a4e 4 + .... In the famous paper |2l| Dyson has argued that this expansion cannot be convergent: 
if this series were convergent, F(e 2 ) would be an analytic function of e at e = 0, and so F(— e 2 ) would 
also have a convergent expansion at sufficiently small e. However, the negative coupling constant e 2 < 
corresponds to the world where the like charges attract each other. The vacuum in this world is unstable 
with respect to the spontaneous creation of very large numbers of particles: the vacuum with electrons in 
one region of space and positrons in another has a lower energy than an "empty" vacuum. Dyson argued 
that the divergence of the power series would manifest itself at high orders ~ 1/e 2 and would physically 
describe the virtual processes of multi-particle production. 

Lipatov [^] has demonstrated explicitly that the degree of divergence of the perturbative expansion 
in scalar theories with dimensionless negative coupling is determined by spherically-symmetrical solutions 
of the corresponding classical equations in Euclidean space. For the same reason, classical solutions also 
dominate large orders of perturbation series if the expansion is performed around an unstable vacuum 123] . 
The action of the classical Euclidean solution determines the probability of tunneling from the unstable to 
the stable vacuum. These methods have been applied to obtain the asymptotic estimates of high orders of 
perturbation theory in scalar electrodynamics [^4] and for a Yang-Mills field coupled to a scalar field 2S]. 



In QCD, we are dealing with the situation in which the perturbative expansion is made around an 
unstable vacuum - the structure of the true QCD vacuum state is not reproduced by perturbation theory. 
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Therefore, we expect that the classical Euclidean solutions of QCD will contribute to the scattering am- 
plitudes, and will describe the processes of multi-particle production. The classical Euclidean solutions - 
instantons [26] - have been found in QCD, as well as in N=4 SYM (see reviews 27, 28, [23] and references 



therein). It has been demonstrated that the instanton "ladders" lead to large multiparticle production 



cross section in QCD [30| at high energy (for instanton- induced effects in high-energy scattering, see also 
[31, 32, 33, 34, 35 1 and references therein). 

All of the instanton-based field-theoretical treatments are restricted to the domain of weak coupling. 
However it appears that instantons also play an important role in supergravity that is related through 
the AdS/CFT correspondence to the strong coupling regime of N=4 SYM. String theory (that reduces 
to supergravity at weak string coupling g s ) possesses solitonic D-branes that are the classical solutions to 
the supergravity equations of motion. If the Euclidean world-volume of a brane is wrapped around an 



appropriate compact manifold, we obtain a so-called D-instanton [36, 37 1 that induces a non-perturbative 
gravitational interaction involving the axion and dilaton fields (see review [38] and the next section). 

The main idea of this paper is to evaluate the multiparticle production that originates from the D- 
instantons in weak gravity. We show that this production is not only large but in the wide region of 
large g s N c generates the total cross section that increases as a function of energy. We argue that the 
main contribution to these processes stems from so-called multiperipheral diagrams of Fig. that we will 
evaluate in N=4 SYM. 

The paper is organized as follows. We introduce all notations and give a general discussion of the 
instanton contributions in weakly coupled supergravity in the next section. In section 3 we calculate the 
diagrams of Fig. Q-type and discuss the energy dependence of the cross section. We show that in spite of 
the fact that the exchange of scalar particles leads to a oc l/s 2 the sum of these diagrams gives a oc s 
with A > 2 in the wide range of strong coupling A. 

2. N=4 SYM in the strong coupling regime: general approach 
2.1 Supergravity (boson sector) 

N=4 SYM at strong coupling is related via the AdS/CFT correspondence Q to the supergravity in ten 
dimensions with the following action (in Euclidean formulation) pSj] 1 : 

_ 5 boson = J d iv x jg- {n - (d^Xd^) + I e ^(d M a)(^a)} (2.1) 
The action of Eq. (|2.1|) describes the interaction of bosonic fields that are given by ten-dimensional metric 



g^ v (gravitons), the dilaton (p and the axion a (that are scalar and pseudoscalar, respectively). In Eq. (2.1) 
1Z is constructed from the Ricci tensor and describes the graviton interactions, while the interaction of the 
dilatons and axions is shown explicitly. 



lr The action of Eq. ( J2.l[ ) is written for ten dimrensional flat space (1Z10). For AdSs the expression is more complicated in 
its gravity part but the instanton contribution remains to be the same (see Ref. [B8|). 
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The first message from Eq. (2.1) is that graviton cannot be produced due to interaction of two dilatons 



(see diagram of Fig. |2|-c which is equal to zero) . The second one is that the large multiparticle production 
can originate only from the dilaton and axion interactions. As we have discussed in the Introduction, this 
multiparticle production can be essential only if the classical solution for the dilatons and axions exists. 

The equation of motion for the action of Eq. (|2.1|) has the following form: 



n pv = i (d^)(d u <P) - ~ e 2 * {d»a){d v a); (e 2 * d^aj = 0; V 2 + e 2 * (d^a)(d»a) = (2.2) 



2.2 AdS 5 X S 5 space 

We start with some basic facts about AdS$ which is defined as the hypersurface embedded in six-dimensional 
flat space-time by the equation (with d=4) 



-Y\ + Y 2 + J2 X i = -L 2 ( 2 -3) 
i=i 

with curvature R = —d(d — 1)L 2 = —12L 2 . Introducing new coordinates 

L Xi L 2 



Y + F-i ' ■ Y + y_! 
we reduce the introduced metric to the following form 



(2.4) 



ds 2 = ^ (dz 2 + = % {dz 2 + dx 2 ) (2.5) 

z V i=i/ 

where x is four dimensional vector. Using Eq. (^1]) we can construct the invariant volume element of AdS^, 
namely 

J d 4 xdz^ = L 5 J d 4 x^r (2.6) 

The metric in the full ten-dimensional AdS§ x S 5 space is 
L 2 

ds 2 = — (dx^dx^ + dyidyi) where [i = 1,...,4 and i = 5, . . . , 10; z 2 = yiyf, (2-7) 
2.3 D-instanton 

The classical (instanton) solution can be obtained from Eq. ( |2.2[) by setting 1Z = 0. The field equations 
then become 

±(d^)(d u 0) = ie 2 ^a)(GU); V M (e 2< ^a) =0; V 2 + e 2< *> ( d tl a)(& 1 a) = (2.8) 



- 5 - 



where all indices are contracted with the ten-dimensional metric of Eq. ( |2.7| ). 
It is easy to see that the solution for the axion field equation is 

a - Ooo = ± (e~* - e - * 00 ) (2.9) 

where both Oqo and (j)^ are the integration constant for the fields at infinity. 
For the dilaton fields we obtain the equation 

g^V^e* = ^dp(yggrdve*) = (2.10) 
Comparing with the equation for the propagator of the scalar massless field in AdS^ [4C] one can see that 



^ = e 0oo + G{x ,z ;x,z) = g s + G(x ,z ;x,z) (2.11) 

where G is a propagator from the point of the location of the instanton (xq, z) to the point of observation 
(x, z). Eq. fl2.9|) gives the instanton solution for the axion field. The general form of G (xq, Zq;x, z) can be 
found in Ref pOfl . Comparing Eq. ( 2.11| ) with the instanton field in the weak coupling regime we see that 



zq plays a role of the instanton size while xq is the position of the instanton. The invariant volume element 
of AdSs (see Eq. (pl3|)) gives us the distribution function n(z) of the instanton with respect to its size 2 , 
n(z) oc l/z 5 . The contribution of the D-instanton to the action is equal to 

2-7T 

•Sb-inst = — = Es P ht = K <y212 ^ 

and can be obtained by direct integration of the field (see Eq. ( |2.11| ) and Ref. p8| ) . The typical time t in 
Eq. ( [2.12|) turns out to be t = z (z is the instanton size). Therefore, the estimate for the typical sphaleron 
energy for the instanton is 

2ir k 

E sp h = = — (2.13) 

9s z z 

3. D-instantons and multiparticle production 

3.1 The general form of the multiperipheral (ladder) diagram 



The multiperipheral diagram of Fig. [T] can be written in more detailed form using Fig. ||. This figure 
shows the ingredients that we need to calculate the amplitude: the correlation function 

(4>inst(xi, Zi) 4>i ns t(xi + i, z^+i )) which is equal to the propagator of the scalar particle in our theory G (a?i+i, Xi, Zi) 
; the vertex for production of the n-particles r (2 — > n) from one instanton; and the probability to find one 
instanton which is equal to 

J ^d'x^Ts = J ^d* Xl e-* (3.1) 

2 Starting from this formula and below we assume that L — 1 to make all expressions simpler. In the final answer we will 
put back the needed power of L. 
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z h 



{3, Zh) 

* 



(21.*)* 



G(x 1 ,z 1 ;32,z 2 ) 



i / 
i / 




The amplitude of the diagram in Fig. ^ 
has the following form in the coordinate rep- 
resentation 

A(x,z h ;x n+ i,z h ) = (3.2) 
dzi 



G (x, z h ; xi,zi) — ^ d 4 Xi e a. 
T{2 



i=l 



A A A\ A AA A A A AAA A A 



unparticles 
Boundary 



Figure 3: The amplitude of multiparticle production in N=4 
SYM in coordinate representation. Dashed lines denote the dila- 
tons and axions, the arrows stand for produced unparticles. 



n i) G (Xi, ZijXi+l, Z; l+ i) 

where the vertex T (2 — > rij) can be calculated 
from the Green's function 

G(x k ,z k \ki,...,k ni ) = (3.3) 

^instixki Zk\x, z) $inst(?Ci) Z%\x , Z ) 

J\ 4>inst{x k o)); 
i=l 



(ri,z—* 0) are the coordinates of the produced particles and (...) denotes the averaging over the position 
of the instanton. The produced particles will be detected at the boundary at large values of afj 3 

The general rules for finding the vertex from such kind of the Green's function is the following (see for 
example Ref. pOp ): go to momentum representation for produced particle, multiply by inverse propagator 
G~ l (ki) , in a such amputated Green's function put particles on the mass shell ( say kf = 0, as in 
the case of gluons). However, we can find the sum over all produced particles (quarks and gluons) in a 
more elegant way by just using the fact that we need to calculate the cross section which is proportional 
to A (x, z n \ x n+ i, Zh) x A* (x, Zh] x n+ i, Zh) and, therefore, what we actually need to know is the Green's 
function (O|0j nst ( 

Xki Zkl^ij Z 

0) 4>i ns t(xk, Zk\ri, z — > 0) 1 0) where Zk, x k and r are coordinates in the complex 
conjugated amplitude. The explicit form for 4>i ns t(xk, 2fc|f», z — ► 0) follows directly from Witten's formula 



</>in8t(Xk,Zk\n,Z -> 0) 



2vr T(4) 



2vr r(4) z\ 



g s ^T(2) ( z 2 + 



Xk 



g s vr 2 r(2)rf 



(3.4) 



We are interested only in the fields at large distances since we consider the scattering amplitude and 
produced particles can be measured only at large distances. Thus the Green's function 
; Zk | r j . z 

0)4>inst(xk — ► 0)|0) can be reduced to ifi\<j>%nst{^i) 4>instiTi)\Q) which can be 
evaluated using the idea of unparticles suggested in Ref. [41]. Namely, 



(0| 4>inst(ri) 4>inst{ri) |0) 



-ik-{f—r) 



\mk)?p{k 2 



d A k 



(3.5) 



(2vr) 4 

where \k) is a state of the unparticle with the four-momentum k^ produced from the vacuum by <j>i ns t- Due 
to the scale invariance 

\(0\4>\k)\ 2 p(k 2 ) = A 2 d=4 Q(k )Q(k 2 )(k 2 ) d - 2 (3.6) 



We want to recall that x denotes the four-vector. 
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Eq. (3.6) sums the emission of quarks and gluons in the most economic way demonstrating one of the 



best features of N=4 SYM approach. Using the dispersion relation one obtains from Eq. ( |3.6| ) that the 
propagator of the unparticle is equal to D unpart ide{k 2 ) = k 4 ln(— k 2 ). We can easily find out the coefficient 
A4 in Eq. (|3.6|) using the following equation 

M = n (VT\ / d 4 re ik - r <f>inst(z;r) = - — (3.7) 

D unp article{k Z ) J Qs 44! 

where z is the instanton size. 



Armed with Eq. (3.5) - Eq. ( |3.7| ) we can explain our pictures in Fig. [l] and Fig. |3|. The dashed line 



in these figures denote the dilaton(axion) fields induced by the instanton. Each of such fields creates on 
the boundary a cluster of produced particles (gluons and quarks); this cluster can be described as an 
unparticle. Using Eq. ( |3.5| ) we can calculate the total inelastic cross section by just summing over the 
number of unparticles without deciphering the content of this cluster in terms of quarks and gluons. 

3.2 Summing over produced unparticles 

From Eq. ( |3.3[ ) one can see that we necoordinateed to sum over produced unparticles since the vertex 
r(2 — » n) oc <j>7 ns t- Eor the production of unparticles with fixed total mass q we need to find 

?-!>] n\]jAl Q(k ^)e(kf)(k 2 ) 2 -^ 

J3.8) 

The factor n! takes into account the possibility that two arbitrary (j)i ns t can enter in (0| (pinsti^i) ( Pinst{^i) |0). 
Replacing 

(2tt) 4 5 (4) [q - k ^J -» J d 4 x exp(i(q-^2ki)-x} (3.9) 

we reduce Eq. ( |3.8| ) to 

S n (g) = A 2n n\ J d 4 xe^ U jj^L ( k 2 ) 2 Q(k ) @(k 2 ) e' its \ (3.10) 

Using the identity D x ex.p(—ik ■ x) = —k 2 ex.p(—ik ■ x) we can rewrite {. . . } in the following way 

{...} = U 2 x l(x) with Z(x) = r^e-bozo [ k \-*W\*\cose dcose \kW 

Jo 27r Jo (2iry 

where \k\ and |x| are the length of the three- vectors k and x repectively. 
After integration over cos# and |fc| we obtain 

T( \ _ f°° dk _ kQX0 k \x\ cos(fc |x|) - sin(fc |a|) 
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Integrating over ko and calculating 0%Z(x) we get 



Ul l(x) 



24-8 



1 



(3.13) 



7T 3 (x2-|x| 2 ) 4 

In all calculations above we assumed that = (q, 0, 0, 0). Using Eq. ( 3.13; ) we have for £ n (<?) the following 
expression 



£n(<?) 



A\ n n\ 



8 • 24 



7T- 



-iA| n n! 



• 24 



7T' 



j4 

axe 



47T 



r 8n 



<£co ■ 



1 vr 3 / 2 r(4n- 3/2) 



,8n-3 



'0 



T(4n) 



A 2 ™ re! 



8 • 24 



7T" 



7r 5 / 2 r(4n-3/2) 8n _ 4 
T(4n + 1) T(8n-3) q 



From Eq. fl3~lj ) we can obtain 

oo oo 



n=l 



n=l 



8-24 



7T 



4 2 

3 A 4 



r(n + l)T(4n-3/2) 
T(4n) r(8n-3) 



!n-4 



(3.14) 



(3.15) 



The sum of Eq. ( pUD 

converges well since at large re the general term decreases as n ' and leads 
to a function that increases with the mass of produced particles q. However, this increase cannot continue 
to large masses due to the unitarity corrections [43] (or due to the quantum effects p4| ). One can see that 
the instanton-induced vertex r(2 — > n) has correction of the type shown in Fig. || and the resulting T,(q) 
cannot be given by Eq. ( 3.15| ) but instead is equal to 



ImT l (q) 



q— >oo 



(1 - 1(g) Re^{q)) 2 + (10 /mS( g )) 2 







(3.16) 



where the imaginary part of Ti(q) is given by Eq. ( |3.15[ ) and (g) stands for all needed integrations; we have 
introduced the notation / for the following integral 



dz 4 ^lE. 
— =- a x e as 

z 5 



dz 



d A x e~ K 



(3.17) 



As will become clear in the next subsection, we actually will need to calculate f dq 2 ImT l (q). The 
integrand here is expected to rapidly fall off above q = E sp h (see Eq. ( 2.13] )) that is the height of the 
sphaleron barrier setting a typical energy scale for the instanton-induced processes. It has been shown 
explicitly |43|, 44 1 that at E > E sp h S(q>) indeed starts to fall off. 

Finally we obtain 

r(re + l)T(4n- 3/2) 

n=l 



T(4re) (4n- l)T(8re-3) 



E 



8n-2 
sph 



(3.18) 
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Figure 4: The unitarity correction to the instanton induced vertex T(2 — > n). 



Using Eq. ( 2.13j ) and the exprerssion for A4 from Eq. ( |3.7| ) we can rewrite Eq. ( 3.18j ) in more convenient 



for for further calculations, namely 



%h rlr?Y(„\ = V r(n + l)r(4n-3/2) / 1 \ n 10re _ 2 

= vr/ ^ ^ fTI^ (4n - 1) f (8 n - 3) I 48^ K (3 ' 19) 



3.3 Multiperipheral kinematics and the interaction between instantons 



\ * r 

ih; gi.±)^ _ J 

\ 1 
\ 1 

\ -x < 

\ 1 
a(oj,A;9i,x)\ / 



\ 

\ 



In this section we recall the calculation of the simple 
multiperipheral ladder diagram of Fig. [| for scalar particles 
in the momentum representation. Since we are interested in 
the total cross section the measure of integration over pro- 
duced particles with momenta k{ = q%—x — qi has the form 
5{k 2 — m 2 ) d 4 fc/(27r) 4 . For momenta qi we use the Sudakov 
decomposition, namely, 



?i(o„,/3„;«„,±)\ / 

\ Y.I 



Figure 5: The simplest multiperipheral dia- ~ a iV\^ + 0iP2,» + %±.;n (3.20) 

gram for scalar particle interaction. Crosses in- 
dicate that the particle is on mass shell. It tumg Qut that the max i ma i contribution for the diagrams 

of Fig. |||-type stems from the following region of integration 

with respect to aj and 

a x > a 2 > . . . » m > . . . > a n ; fix <C <C . . . < Pi < . . . <C /3„; (3.21) 

One can see that the integration can be performed using momenta % with the phase space doti dfiis d 2 qi j ±/(2ir) 4 
where s = (px + P2) 2 and 

q 2 = aipis + q 2 ± = q 2 ± + O (m 2 ^-) ; (3.22) 

,2 _ /„ . o , „. a \ „ iu2 _ „. a „ ,2 , , n /\2 «i 



fcj = + A) s — kf± = a-ifii+xs - k i± + 0[m 



Cii+X 
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Therefore, the diagram of Fig. |5| leads to the following expression 
A 2 n (Fig.§ = 



f d 



4 TT I doii d Qh-L 9 2 x I 2 2 \ x ( a 2 2 \ 

9 11 J ( q 2 +m2)2 5 5{a lS -m - g 1>± ) 5 (/3 n * - m - g n J 



1 g ( f dai d ' q% ± don 9 



n\ s \J an=m 2 /s a>i (2vr) 3 a { (gf ± + m 2 ) 2 ^ 
5 4 ln n (s/m 2 ) / d 2 g ij± 5 2 V 



(2vr) 3 {q 2 ± + m 2 ^ 2 



(3.23) 



Summing over n we get the total cross section: 



1 00 A— 2 

= - E ^ = 5 4 (4) " . ( 3 - 24 ) 



n=l 



where 



(3.25) 



(2tt) 3 (g 2 ± + m 2 ) 2 " 

The expression for A can be re-written in the impact parameter (6) representation: replacing the particle 
propagator in q -representation (G(q) = l/(q± + m 2 )) by the propagator 



we obtain that 



The diagram of Fig. || is more complicated. First, we do not produce one particle but the bunch of particles 
that we represent by an unparticle. However using the dispersion relation for £(<?) we can replace 

1 1 f ZmS(M) dM 2 



rn 



2 -k 2 2tt 



J Im ^[ M ^ 2 M ; 5{m 2 -k 2 ) — > I dM 2 5{M 2 -k 2 ) ImZ(M) (3.28) 



Integrating over /3j we obtain an extra factor J* dM 2 JroE(M). 

A natural question arises about what should we use instead of G(b) in Eq. ( |3.27| ). As has been discussed 
above, the instanton field is equal to pS[ | 

27T 

<t>inst{zk,x k ;ZiX) = — G(z k ,x k ;z,x) (3.29) 
5s 
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where G (z k ,x k ; z, x) is the dilaton propagator in AdS§ x S 5 theory. Therefore to find the Green's function 
for the interaction of two instantons we need to evaluate the following product: 



<f>i ns t(z k ,x k ;z,x)G 1 (z k ,x k ;z,x) x G (z, x; z' , x) x G 1 (V, x'; z k _ x , x k _i) <f>i ns t(zk-i, x k -ii z',x') 



2ttV 



J G(z h ,x k ;z k -x,x k -i) = K 2 G(z fc ,a; fc ;z A .„ 1 ,x fc _i) 
9s J 



(3.30) 



It is clear that this procedure is just the same as we did in Eq. ( |3.5| )-Eq. fl3.7| ) for unparticles that have 
been produced on the boundary. 

As we have learned the most convenient way to factor out the energy dependence is to perform all our 
calculation in the mixed representation: two energy variables (on and f3i (or qi )+ and qi~) and coordinates. 



In AdS§ such propagator has been found (see Refs.p, 11, |i4[): 



G(z /c ,x fc ;2; fc _i,a; fc _i) = G 3 (u) = — 



1 



+ u + vXit + 2)} V u ( u + 2 ) 



with 



(3.31) 



(z k - z fc _i) 2 + b 2 



(3.32) 



2 z k z k -\ 

where b = x k -i t ± — x kj ± is the impact parameter which is a variable conjugated to qi in Fig. [5|. 

For the vertex T(2 — > n) in Fig. ^ we need to collect factors from Eq. ( 3.3C| ) and Eq. ( |3.17| ). Doing so 



we obtain that we need to replace g in Eq. (3.23) by 



dz k 



(3.33) 



Concluding this section we can claim that \A\ 2 for the amplitude of Eq. (|3.2| ) can be written in the 
multiperipheral form of Eq. ( 3.22j ) and Eq. ( |3.23| ) with the following short glossary 



5{m 2 - k 2 ) 
G{b) (see Eq. flpf)) 
9 



dq 2 E(g) (see Eq. ( p9p ) \ 
G{u) (see Eq. ( pl|) ) ; 

^ (see Eq. O) 



(3.34) 
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4. The energy dependence of the multiparticle production cross section 
4.1 Equation and its solution 



z') 



For the imaginary part of the scattering amplitude we can write 
down the equation which is illustrated in Fig. |6[ It has a simple form, 
namely 



d$(z,z) /2tt\ 4 _ 2 22l 
— — e ss 



dins 



\9s J 



(4.1) 



where the kernel K is given by 



K(z,z;z',z') = d 2 bG(z,z';b) G(z,z';b) dq 2 S(g) (4.2) 



Figure 6: The equation for the scat- As we have geen in Eq (g3g[) J ^2 s ^ ig pr . opor . tiona l to z 'z'. In- 

tering amplitude. The black crosses troducing J dq 2 s(?) = ^ ^( K ) we can rewrite Eq. Q in the form 
denote the instanton vertices. is 

given by Eq. (|3T5|) while propagators d$(z,z) 4 _ 2k ~ f dz'dz' 

G are written in Eq. (|3.3l|) . rflns - K e E(k) 



z' 4 z' 4 



K(z,z;z',z) $(z',z) (4.3) 



with 



d 2 bG(z,z';b) G(z,z r ;b) 



(4.4) 



The general form of this kernel is rather complicated but in the region where z ^> z' and z 3> z' the 
integral over b can be taken explicitly and it leads to 

(45) 

where 

z\ -1 + 8x 2 - 8x 6 + x 8 + 12x 4 ln(x 2N 



Fix = - = ■ — r — (4.6) 

V z) 4(x 2 - l) 5 V ; 

with t = z z and t' = z' z' . Therefore, to study the main features of the kernel K of Eq. ( [4.4| ) we can use 
the simplified formula 

K [z, z; z', z') = i- | £ F(x) 6 (t - t') + £ F(x') 6 (f - t) } (4.7) 

where x' = z 1 jz! . 
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G(z h ,z';b) 



We start the analysis of the solution to Eq. flbl| ) by calculating 
the first iteration of this equation given by Fig. [/]. For the region of 
integration over t' < we have the contribution 



Figure 7: The simplest diagram of 
the instanton contribution (the first 



A 2 (Fig. 0; t h , 
Ms/go) 

2vrs 



/*h r z h i +6 i/6 

dz'y d - z '-F 2 ( Xh )^- 



4 - 

k e 



v 7 3 v 



(4.8 



interation of Eq. (4.1)). The black 
crosses denote the instanton vertices. 



In the case when t' > th one can see that 



E(g) is given by Eq. (3.15) while prop- 



agators G are written in Eq. (3.31) 



A 2 (Fig.%t h ,x h ) 
ln(s/s ) 



2-ks 

Mf/fo) 
2tts 



k e 



E(«) 



1 / 6 /' 6 

z h J Zh L l h 



dz 



k e 



- 2 «£( K )±F 2 (x h )±- 



(4.9) 



The lesson from these simple estimates is clear: the typical z' and z' are about Zh and z^, respectively. 
Actually, we need to integrate over Zh and Zh separately for target and projectile and for high energy 
scattering Zh — Zh (see Refs.|^, 0). However for our purposes the simple estimates given above suffice. 

Comparing the sum of Eq. (|4.8|) and Eq. (|4.9| ) with the diagram without the instanton which in our notation 
is 



A 2 (without instanton) = ty L F(x)/s 

we can conclude that 

A 2 (Fig. ^) = A 2 (without instanton) ln(s/«o) I /; ' 
= A 2 (without instanton) A ln(s/«o) 

with 



(4.10) 



2k 



E(«) 



10 

21 ^ >/ 



A 



2vr V ; 21 t 2 



(4.11) 



(4.12) 



In this estimate we assumed that the final answer will be A 2 oc s A_1 . In spite of the fact that these simple 
estimates are very instructive we should not expect them to give a reliable quantitative result since we 
made a strong assumption that z and z in the target and projectile are the same. 

For the solution of Eq. (|4.1|) we assume that 



/ 



dr$(z, z) 



z 2 z 2 $(zz\s) 



(4.13) 
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12 



Figure 8: A(/t) versus k, for = 2 a' 2 and N c = 3. 



where r = z/z, and looking for 

For cj) (t) Eq. ( |4.1| ) has the form 

2 d(j)(t;s) 



<&(t = zz) = 4>(t h = z h z h ) s A 



dins 



where 



K{t) 



I e~'^ / dt'K(t') <£(*'; s) 
9s J 



Eq. fl4.15| ) stems from Eq. ( |4.l| ) after integration over r of both parts of this equation. 



(4.14) 

(4.15) 
(4.16) 



The integrals over z' and z' can be taken by switching to a new variable y = 2zz' yu + 1 + \/u(u + 2) 
With this variable the integrals can be easily evaluated: 



dy 

4^ 7 2( ^ +62) y 2 (y-2z 2 ) ' 4vr z 2 (z 2 + 6 2 ) 



(4.17) 



and 



= / d r [ d 2 bl(z;b) I(z;b) = — f dr , 1 In (1/r) = — — = — 
w y y v , ; v > ; 167r y r (l-r 2 ) v 7 ; 16vr 12 1612 



Finally, we obtain for A the following equation 



^ J_ ... i „-2k 

16 12 t 2 



^K 4 e^ K S( K ) 
/l 



(4.18) 



(4.19) 
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4.2 The Pomeron intercept 

Eq. (|4.15|) gives us the energy dependence of the cross section, namely, 

(^inelastic OC S A_2 (4.20) 

Therefore our mechanism of multiparticle production can be important at high energies only if A exceeds 
2. At first sight A has a dimension since it is proportional to 1/zf. This is the artifact of our simplification 
since we put L = 1. Restoring L explicitly, we get the following expression for A: 

a = f « 4 ^m = i^f" 4 ^"" = £ n 4 k3 ^ [k) (421) 

To take account of the instanton interactions with axions, we need to multiply Eq. ( [4.21| ) by factor 
of 4 ( see Eq. ( |2.1|) ). Indeed, we have four transitions in the t-channel: 2d — ► 2d, 2d — > 2a, 2a — ► 2a and 
2a — > 2d (d and a denote the dilaton and axion respectively). 

In Eq. ( f4.21| ) we used that A = g\ M N c = ocym 4vr N c = g s 4ttN c = 8tt 2 N c /k and L 4 = a' 2 A, where 



a' is the slope of the Reggeon trajectory, g s is a string coupling, and gyMdinstmp.tar is the YM coupling. 

In Fig. H we plot A(/t) assuming that z\ = 0.5 a' , and N c = 3. We believe that this case can 
resemble the real world (see discussion below). One can see that for k > 4 A(/t) > 2 and therefore the 
instanton-induced multiparticle production has a cross section that increases with energy. 

At k = 2 A = 120 but the value of g s ~ 3 is rather large for the approach that can be reduced to weak 
gravity. We recall that this approach corresponds g s < 1 but g s N c > 1. From Fig. || one can see that at 
k = 8 we have g s ~ 0.75 while g s N c = 2.25 > 1. For this value of g s A ~ 30 and the graviton interaction 
generates (2/VX) w 15 - 20% of the cross section. Our approach can be trusted if g s < 1 but g s -/V c > 1; 
this translates into the window for k in Fig. ^ is 6 < k < 18. One can see that for values 6 < k < 8 
which are in this window A — 2 > leading to the cross section for multiparticle production that increases 
with energy. 

These estimates show that we can obtain phenomenologically reasonable values of A. However, these 
values crucially depend on the value for Zh which gives us a new scale for high energy interaction. We 
cannot calculate the wave functions of the observed hadrons in N=4 SYM (remember that this theory 
does not possess confinement). However in our estimates we implicitly assumed that the hadron wave 
functions are peaked at z = Zh- As has been seen from our calculation of the Pomeron intercept and as 
we will demonstrate in the next section calculating the slope of the new Pomeron trajectory, z^ is our key 
phenomenological parameter that determines all features of the Pomeron. We choose z\ = 0.5 ol since 
it gives reasonable Pomeron intercept and a' p ~ 0.5 a' as it should be for the closed string in the string 
theory that corresponds to N=4 SYM in the limit A — > oo. 

4.3 The Pomeron slope 



Eq. (44J allows us to find the slope of "Pomeron" trajectory A(t) = A + a' P t where t = —Q is the 



momentum transferred along the Pomeron. However we need to rewrite the equation for t ^ 0; this leads 
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to 



with 



d®(z,z) 
din s 



2ir\ f dz' dz' / , „ . T , , ,-. 

e 2 »- / K(z,z;z',z'\Q) <S>(z',z>) 



9s J 



z J z 



K(z,z;z',z'\Q) = d 2 b G(z,z';b) G(z,z';b) dq 2 E(g) 



(4.22) 



(4.23) 



Actually for calculation of the slope we need only to find the second term in the expansion of K (z, z; z', z'\Q) 
with respect to Q, namely, 

K (z,z;z',z'\Q) = (4.24) 
= J d 2 b G (z,z';b) G(z,z';b) J dq 2 E(g) + ^Q 2 J b 2 d 2 b G (z, z';b) G (z,z';b) j ' dq 2 S(g) 
= K\ {z, z ; z', z') + Q 2 K 2 (z,z;z',z) 



K\ is the kernel of Eq. (44) for which we have the explicit form given by Eq. ( |4.5| ) for z 3> z' and z 3> z'. 
In the same kinematic region K 2 has the form 



K 2 (t) = dr j b 2 d 2 bl(z;b) I(z;b) 

Using Eq. ( flgg) we obtain the expression for the slope of the Pomeron trajectory: 



(4.25) 



a' P = (K 2 (t h ) fK x (t h ) ) A = 0.23 t h A (4.26) 

The numerical coefficient is the result of the numerical evaluation of the integral of Eq. ( |4.25| ) . Using our 
value for = 0.5a' and A rs 4 which we have for N c = 3, we obtain a' p « 0.5 GeV~ 2 . 

4.4 Corrections 



a^(SS, ft) 



9 , ,92-9 



_\ -x f _ J 



93(03, ft; 93,±) \ 



9i(a„,/3„;g^±) \ / 

m \ y/ 



Figure 9: The correction to the amplitude due to Figure 10: The correction to the amplitude in the 
exchange of two dilatons. All notation as in Fig. ^| simple multiperipheral diagram. 

In this subsection we will discuss the possible corrections to the simple multiperipheral diagram of 
Fig. ||. The main corrections stem from the two-dilaton exchange (see Fig. |9|). It easy to see that these 
corrections change the value of A. 
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Indeed, let us consider the simple multiperipheral diagram of Fig. 10. The extra dilaton in t channel 
will generate the imaginary amplitude and therefore only the diagram with the exchange of two extra 
dilatons can contribute to the total cross section. To evaluate it we need to perform the extra integration 
over d 4 q = da q d(3 q s d 2 q± where 

% = a q Pi,n + (3 q P2,n + q± lf i (4.27) 

We close the contour of integration over a q and (3 q over poles: (qi — q) 2 = and (q — q 3 ) 2 = 0. In doing 
so we obtain 



/ 



sda a dB a d 2 q\ 1 1 1 , m? , „ m? 

\ q t S2 7 \2 = 27r i —5— with a i = - 7^- and & = — 428 

2m (q 1 - q) 2 (q - q 3 ) 2 a x /3 3 s f3 3 s «is 



where m 2 = m 2 + (qi t ± — 92,±) 2 or m 2 = m 2 + (q 3j ± — q2,±) 2 - ^From 5 ((qi — (ft) 2 ) and 5 (((72 — Q3) 2 ) we 
know (see section 3.3) that we have 

f 2 - f with a x = ^- (4.29) 



P 1 /?2 J fa /?3 fts 



Multiplying this integral by Eq. ( 4.28 ) one can see that we have 



M / f j f 1 K / dft / § (430) 



Including the exchange of the extra dilaton one can see that the integral has the following structure 

dfo f A 



W / fodfo / — » / ^ (4.31) 



Therefore, instead of two logarithms from the integration over (3% and (5 3 (see Eq. ( 4.29| )) we obtain only 
one logarithm. 



To take these corrections into account we need to change the kernel in Eq. (4.1), 

K (z,z;z',z') — ► K (z, z; z', z') 

+ (J)' e~ 2 % J dq 2 J ^ |J G 2 (z, z"; b) G 2 (z, z"; b) d 2 b K (z" , z"; z>, z>) (4.32) 

The numerical estimates show that the second term can reach about 1 -f- 5 % of the first one. This 
smallness stems from the ratio J ^pr J G 4 d 2 b/ J G 2 d 2 b ~ 3 10 -4 [x" = z" /z"). In these estimates we used 
the simplified version of the kernel given by Eq. 

5. Conclusions 

In this paper we have demonstrated that the D-instanton induced processes can be a dominant source of 
multiparticle production in the range of the string coupling constant g s < 1 but g s N c > 1. The cross section 
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of the instanton-induced multiparticle production increases with energy and the resulting picture of high 
energy scattering in strongly coupled N=4 SYM looks as follows. The weakly coupled graviton generates 
the elastic amplitude and the corresponding part of the total cross section, while the D-instanton induced 
interactions of dilatons and axions are reponsible for the mutiparticle production processes. This two 
component picture leads to the structure of the high energy interaction which is close to the experimental 
observations as well as to the QCD predictions. 

We would like to stress that we are forced to look at the dilaton and axion sector for the source of 
multiparticle production: as we have discussed above, the AdS/CFT correspondence reduces the N=4 
SYM in the strong coupling limit to the action of Eq. ( |2.l[) in which the gravitons decouple from dilatons 
and axions. Since the graviton interactions in weakly coupled gravity can only induce elastic interactions 
(responsible for the Newton's law), the only possibility to describe the multi-particle production is to 
invoke the dilaton and axion sector. We hope that the proposed here D-instanton-induced mechanism of 
multiparticle production in strongly coupled N=4 SYM provides a resolution of the following dilemma: 
the multiparticle final state (the strongly coupled plasma) successfully described by this theory seemingly 
could not be produced in the scattering processes dominated by the elastic graviton exchange. 

Our estimates of the value of the intercept of the new 
inelastic Pomeron should be taken with caution: the estimate 
for z\ is not reliable and the equation was solved with the 



0.35 



0.25 




approximate kernel of Eq. (4.7) instead of one of the full one 
given by Eq. ( |4.3| ). Our numerical calculations show that 
using the full kernel we obtain A approximately two times 
larger than with the simplified kernel. In Fig. [ll] we plot 
the dependence the maximal and minimal values of 2/V\ 
versus N c that satisfy the conditions :<? s < l,g s N c > 1. For 
these 2 / y/\ the value of A is larger than 2 and the instanton 
induced multiparticle production increases with energy. 

It is interesting to compare the pictures of multiparticle 
production based on QCD instanton ladders at weak cou- 
pling |3(| and on D-instanton ladders in strongly coupled 
N=4 SYM. In the D-instanton picture, we are dealing with 
the production of unparticles (bunches of quanta with con- 
tinuous invariant mass distribution). This is analogous to the production of continuous spectrum of gluons 
from the decay of the QCD sphaleron. Another feature of the D-instanton picture is the appearance of 
a new non-perturbatives scale: z%. It plays the role of the typical size of the instanton in QCD vacuum 
that is taken from lattice simulations. It appears that the typical size of QCD instantons and the size of 

V0.5 a 1 



Figure 11: The maximal and minimum 2/%/A 
versus N c for which A > 2 and g s < 1 while 
9s N c > 1. 



D-instantons in our present estimates are quite close to each other: we have used z^ 



0.2 fm. 



D-instanton can be viewed in the string frame [36| as an Einstein- Rosen wormhole connecting two 
asymptotically Euclidean regions of space-time, with the Ramond-Ramond (RR) charge flowing down the 
throat of the wormhole. It describes a violation of the conservation of a global charge in physical processes 
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[ 36 1 . Since in the AdS/CFT dictionary the RR scalar (or the axion, in our notations throughout this paper) 
of supergravity is dual to the 9 angle of N=4 SYM field theory [45], this non-conservation of RR charge 
describes the change of chirality 4 . D-instantons thus play a role that is very similar to the role played by 
the QCD instantons that violate chirality conservation. 

It has been predicted Ef| that this violation of chirality in hot QCD matter would induce the charge 
asymmetry in high energy heavy ion collisions with respect to the reaction plane (i.e. the quark-gluon 
plasma would develop an electric dipole moment fluctuating on the event-by-event basis). Detailed theory 
of this "chiral magnetic effect" (chiral charge generates an electric current along the direction of external 



magnetic field) has been developed in Refs. [47, 48, [4{j, 50, 51 1. The evidence for the chiral magnetic 
effect from the lattice QCD has been found in [52 
been explored by holographic methods in (5^ 



The behavior of the effect at strong coupling has 
, 55 1 . Very recently, STAR Collaboration has presented 
an observation of charge asymmetries in heavy ion collisions at RHIC [56, 57 1. Our D-instanton-based 
approach suggests an important role for topological effects in high energy collisions. 
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